Phase transition and properties of compact star 
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We investigate the phase transition to a deconfined phase and the consequences in the formation 
of neutron stars. We use the recently proposed effective field theory motivated relativistic mean 
field theory for hadron and the MIT Bag model and color-flavor locked (CFL) phase for the quark 
matter in order to get the appropriate equation of state. The properties of star are then calculated. 
The differences between unpaired and CFL quark matter are discussed. 
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I. INTRODUCTION 



One of the fundamental problem of nuclear physics is 
to understand the behavior of nuclear matter at extreme 
conditions. The study of neutron star provide an im- 
portant information in this regard. Neutron star are 
extremely dense objects. They formed after the gravi- 
tational core collapse of a massive star The envi- 
ronment in central region resembles the early universe, 
except that the temperature (T) is lower. Due to the 
lower temperature and high density, neutron stars are 
presumably unique astrophysical laboratories for a broad 
range of physical phenomena [l|, The composition 
and other properties of neutron stars are depend upon 
the appropriate equation of state (EOS) that describe 
its crust and interior region Q. The crust part of neu- 
tron star, where the density is comparable to saturation 
density (po) of symmetric nuclear matter, is adequately 
described by hadronic matter. But the interior region 
where density p is of the order of 5 to 10 times of po, a 
phase transition occurs from hadronic to quark matter. 
But this is not well understood, whether the interface re- 
gion, it is a quark and/or mixed matter. It is reasonable 
to assume that the central part is dense enough, so that 
it can be treated as a pure quark phase. However, the 
boundary part between the quark and hadronic matter 
should be an admixture of quark and baryonic degrees 
of freedom. In the present study, we are interested in 
building the EOS for mixed matter of quark and hadron 
phases. We employ an effective field theory motivated (E- 
RMF) Lagrangian approach including hyperons in order 
to describe the hadron phase. For the quark phase we 
have chosen to use both unpaired quark matter (UQM) 
described by the MIT bag model @, [E H and paired 
quarks described by the color-flavor locked (CFL) phase. 
Recently many authors [1 S El El El El 

141 have dis- 
cussed the possibility that the quark matter is in a color 
superconducting phase, in which quarks near the Fermi 
surface are paired, forming Cooper pairs which condense 
and break the color gauge symmetry [15| . At sufficiently 
high density the favored phase is called CFL, in which 
quarks of all three colors and all three flavors are allowed 
to pair. 

The paper is organized as follows: in section II the E- 



phase for quark matter and mixed phase are described. 
Here, we are trying to build an equation of state for mixed 
matter of hadron and quark phases. We use an E-RMF 
model including hyperons with incompressibility K = 
300 MeV and M*/M = 0.7 of nuclear matter for hadron 
phase. We use Gibbs criteria and chemical equilibrium 
conditions, to built a mixed phase EOS, and then we 
calculate and discuss the properties of star at T=0. The 
calculated results are discussed in section III. In the last 
section, the conclusions are drawn. 



II. FORMALISM 
A. Hadron matter 

In principle, one should use quantum chromodyan- 
mics (QCD), the fundamental theory of strong interac- 
tion, for the complete description of EOS. But it can- 
not be use to describe hadronic matter due to its non- 
perturbative properties. A major breakthrough occurred 
when the concept of effective field theory (EFT) was in- 
troduced and applied to low energy QCD [16j]. The EFT 
for strong interaction at low energy is known as quan- 
tum hadrodynamics (QHD) [H El, El M, S3l- The 
mean field treatment of QHD has been used extensively 
in order to describe properties of nuclear matter [171 |22I | 
and finite nuclei [11 SI [H H]. The degrees of free- 
dom in this theory are nucleons interacting through the 
exchange of iso-scalar scalar a, iso-scalar vector ui, iso- 
vector-vector p and the pseudoscalar n mesons. The nu- 
cleons are considered as Dirac particle moving in classi- 
cal meson fields. The contribution of tt meson is zero 
at mean field level, due to pseudo-spin nature. The 
chiral effective Lagrangian (E-RMF) proposed by Furn- 
stahl, Serot and Tang [2?], [H, [2t| is the extension of 
the standard relativistic mean field (RMF) theory with 
the addition of non-linear scalar- vector and vector- vector 
self interaction. This Lagrangian includes all the non- 
renormalizable couplings consistent with the underly- 
ing symmetries of QCD. Applying the naive dimensional 
analysis [3d |31| and the concept of naturalness one can 
expand the nonlinear Lagrangian and organize it in in- 
creasing powers of the fields and their derivatives and 
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tice, to get a reasonable result, one needs the Lagrangian 
up to 4th order of interaction. In the interior of neu- 
tron star, where the density is very high, other hadronic 
states are produced [H, Ha, [H, Thus, the consid- 



ered model involves the full octet of baryons interacting 
through mesons. The truncated Lagrangian which in- 
cludes the terms up to the fourth order is given by 
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The subscript B = n,p, A, £ and S, denotes for baryons. 
The terms in eqn. (1) with the subscript "B" should be 
interpreted as sum over the states of all baryonic octets. 
The covariant derivative is defined as 



+ iguiB^^ + ig<j>B<t>n + ig P BhBT a p a ^, (2) 



whereas R^ u , and Cl^ are the field tensors 



(3) 



a, 



(4) 



m p are the 



where m B denotes the baryon and m a , m u , 
masses assigned to the meson fields. Using this La- 
grangian, we derive the equation of motion and solved it 
in the mean field approximation self consistently. Here, 
the meson fields are replaced by their classical expecta- 
tion values. The field equations for a, uj and p-meson are 
given by 
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For a baryon species, the scalar density, pss, an d baryon 
density (/9 b) are 
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where 



\/k 2 + m* B is the effective energy and J B 
and /3s are the spin and isospin projection of baryon J5, 
the quantity k B is the Fermi momentum for the baryon, 
to* = m B ~ g aB u is the effective mass, which is solve self- 
consistcntly from equation (5). After obtaining the self- 
consistent fields, the pressure P and total energy density 
e for a given baryon density are 
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As explained earlier, the terms in equations (10) and (11) 
with the subscript "2?" should be interpreted as sum over 
the states of all baryonic octets and 7 = 2 is the spin de- 
generacy. In the above, Pi and Si are lepton pressure 
and energy density respectively, explained in the follow- 
ing subsections. 

For stars in which the strongly interacting particles 
are baryons, the composition is determined by the re- 
quirements of charge neutrality and /3-equilibrium con- 
ditions under the weak processes B\ — > B2 + I + vi and 
B 2 + 1 — > Bi + v\ . After deleptonization, the charge neu- 
trality condition yields 

ftot = $>b(2J b + l)k 3 B /(67r 2 ) + mk?/(3ii 2 ) = , 

B 
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where qs corresponds to the electric charge of baryon 
species B and qi corresponds to the electric charge of 
lepton species I. Since the time scale of a star is ef- 
fectively infinite compared to the weak interaction time 
scale, weak interaction violates strangeness conservation. 
The strangeness quantum number is therefore not con- 
served in a star and the net strangeness is determined by 
the condition of /3-equilibrium which for baryon B is then 
given by hb — fteMn — toMej where [Ib is the chemical 
potential of baryon B and bs its baryon number. Thus 
the chemical potential of any baryon can be obtained 
from the two independent chemical potentials fi n and fi e 
of neutron and electron respectively. 

The lepton Fermi momenta are the positive real solu- 
tions of (fcg + ml) 1 / 2 = fi e and (k 2 + m 2 ) 1 / 2 = (i^ = fj, e . 
The equilibrium composition of the star is obtained by 
solving the set of Eqs. (5)- (7) in conjunction with the 
charge neutrality condition (|12[) at a given total baryonic 
density p = J2 b (2Jb + 1)&b/(67T 2 ); the baryon effective 
masses are obtained self-consistently. 



B. Unpaired Quark Matter 



which the u,d and s quark are degrees of freedom with 
electron. In this model the masses of u and d are set to 
5.0 MeV and strange quark mass is taken to be 150 MeV. 
The chemical equilibrium is given by 



Pd = fJ-s — M« + Me 



(13) 



The fi n and fi e are the two independent chemical poten- 
tials and rest can be written in terms of them as follows: 



Mu = g Mr, 



Md = gM« + 7j Me 
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(14) 

(15) 
(16) 



The pressure for quark flavor /, with /= u,d or s is 
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where kf = (p 2 — to 2 ) 1 / 2 is the fermi momentum. 
The leptons pressure is 



p A dp 



(18) 
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The total pressure is given by 

P = Pi+P q -B (19) 
where B is the bag constant. 



In the central part of neutron star the density is ex- 
pected to high enough that hadronic matter undergoes to 
quark degrees of freedom. In quark phase we employ MIT 
bag model to describe unpaired quark matter [38, 39, 40]. 
The bag model provides a useful phenomenological de- 
scription of quarks being confined inside the hadrons. 
Confinement results from the balance of the bag pres- 
sure on the bag walls from the outside and the pressure 
resulting from the kinetic energy of the quarks inside the 



C. Color Flavor Locked Quark Matter 

In this section, we consider the quark matter phase 
as color flavor locked (CFL) quark paired, in which the 
quark near the Fermi surface form Cooper pairs which 
condense, breaking the color gauge symmetry. In CFL 
phase, all the three colors and three flavors are allowed 
to pair and this is the favored phase at sufficiently high 
density. We describe the CFL phase by using the ther- 
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and the Goldstone boson contribution due to chiral sym- 
metry breaking in the CFL phase is given by 
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and the electron contribution f2/ is given by 
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The quark number densities and the electrical charge 
density carried by the pion is given by 



Pu= Pd= Ps 



2A 2 fi q 



(27) 



where A is gap parameter and its value is 100 MeV, which 
is the typical value considered in the literature. 

The electric charge density carried by the pion conden- 
sate is given by 



{CFL 



,T/M Me 



(28) 



In the above thermodynamic potential, we have neglected 
the contribution due to the kaon condensation which is 
an effect of order m\ and thereby small compared to 
theA 2 /^ 2 contribution to the thermodynamic potential 
for A ~ 100 MeV. 



D. Mixed Phase and Star Properties 

The mixed phase is obtained by applying charge neu- 
trality condition and Gibbs criteria for hadron and quark 
phase. The charge neutrality condition is: 



where p® p and p^ p are the charge density of quark and 
hadron phase, x an d (1 — x) are the volume fraction 
occupied by quark and hadron phase respectively. The 
phase boundary of the coexistence region between the 
hadron phase and quark phase is determined by Gibbs 
criteria. The critical pressure, critical neutron and elec- 
tron chemical potentials are determined by the condi- 
tions: 



PHP,i = PQP,i = Pi,i = n,e 
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(30) 



(31) 



(32) 



These are the chemical, thermal and mechanical equi- 
librium, respectively. The energy and the total baryon 
densities in the mixed phase are: 



e = xe QP + (l-x)s HP +e t , 



and 



XP 



qp + (i-x)p hp - 



(33) 



(34) 



Once the above quantities are obtained, we can construct 
the EOS for mixed phase and consequently compute the 
properties of neutron star. To evaluate the star structure 
we use the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions found in Ref. |42| . They are 



dP 

dr 



G [e + P][M + Airr 3 P] 
7 r — 2GM : 
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where G and M{r) are the gravitational constant and 
enclosed gravitational mass and c = l. For a given EOS, 
these equations can be integrated from the origin as an 
initial value problem for a given choice of the central 
density, (sq). The value of r(= R) at which the pressure 
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III. RESULTS AND DISCUSSIONS 

We obtain the EOS for hadronic matter by changing in- 
compressibility K=300 MeV and effective mass (M*/M) 
to 0.7 of G2 parameter set [43| , It is to be noted that we 
are not getting any mixed phase with original G2 set. So 
we change the incompressibility K from 215 MeV to 300 



MeV and effective mass (M*/M) from 0.664 to 0.7. The 
resultant parameter set satisfies all the nuclear saturation 
properties so that our extrapolation to higher density re- 
mains meaningful. After the above re-adjustment, the 
modified Lagrangian parameters are given in Table I. To 
compare the changes, the original G2 set is also displayed 
in the Table. 



TABLE I: G2 and modified G2 parameter set (G2*) 



Set 


m a /M 


m w /M 


m p /M 


9 s /4tt 


3^/4tt 




K3 


K4 


Co 


Vi 


V2 




G2* 
G2 


0.554 
0.554 


0.833 
0.833 


0.820 
0.820 


0.751 
0.835 


0.936 
1.016 


0.614 
0.755 


1.202 
3.247 


14.981 
0.632 


2.642 
2.642 


0.650 
0.650 


0.110 
0.110 


0.390 
0.390 



We assume that all the hyperons in the octet have the 
same couplings. They are expressed as a ratio to the nu- 

cleon coupling x a = XHa/xNa = \J\, x u = x H u/xnu = 

and x p = XHp/xpjp = yf- Using the modified 
values of G2 set of Table I, the EOS for quark mat- 
ter obtained with unpaired quark model and color flavor 
locked phase for different values of the Bag pressure B. 
In our calculation, we took B* = 170 MeV, 188 MeV for 
UQM and CFL models respectively, and then calculate 
the equation of state. For the CFL, it is to be noted that 
we do not find any mixed phase at = 170 MeV and 
180 MeV. In both the cases we have m u = = 5 MeV 
and m s = 150 MeV. 

In Fig. 1 we plot the resulting EOS for both E- 
RMF+UQM and E-RMF+CFL cases. From Fig. 1 it 
is evident that the mixed phase starts earlier in case of 
UQM whereas CFL predicts the mixed phase at higher 
density approximately at around e = 2 fm -4 . The inclu- 
sion of hyperons softens the EOS [3| of the charge neu- 
tral dense matter. This is clearly seen from the change 
in the slope of Fig. 1 for energy density greater than 
e w 2.2 fm" 4 . 

In Fig. 2, the fractional particle densities, Pi/po for 
baryons, leptons and quarks in E-RMF+UQM are shown. 
From the Fig. 2, we notice the onset of quarks at around 
1.3po and immediately in the vicinity of 4po they are the 
most abundant particle in the matter. With increasing 
density the model predicts pure quark matter after 6p Q . 
At this density and thereafter all the three quarks under 
consideration contributes almost equally to the matter 
density. As evident from the figure that the leptonic con- 
tribution ceases after 4po, which are primarily used up 
to maintain charge neutrality of the matter. The nucle- 
ons constitutes a sizeable population in the matter but 
after they decrease abruptly at around 5po, after which 
its a pure quark phase. Further it can be seen that A 



ticle fraction decreases, finally ends at around 6po- The 
appearance of A is decided purely by the neutron chemi- 
cal potential since it is isospin independent and hence the 
density of A decreases with decrease in neutron density. 
However the contribution of A is quite negligible overall. 

In Fig. 3, the particle population for the baryons, lep- 
tons and quarks are shown as a function of baryon density 
up to 8p in case of E-RMF+CFL for B* = 188MeV. 
From the Fig., it can be seen that the appearance of 
quark starts ~2.5po- It is to be noted that the quarks 
are of equal densities in CFL. The deleptonization in the 
matter occurs at ^5.5po- Similarly, the nucleons follows 
similar trend and at ~5.5po the matter is in pure quark 
phase. Like Fig. 2, A appears at Ap and decrease with 
the nucleons particle fractions and finally ends at 5.5po- 
In CFL case, E~ appears at 1.89 po and ends at 3.5po- 
In both the cases i.e., E-RMF+UQM and E-RMF+CFL 
the difference is due to the different charge neutrality 
conditions. 

The strangeness content in case of high density mat- 
ter is an exciting possibility, which can give important 
insights into some of the most fundamental problems of 
astrophysics and high density behavior. Figure 4 dis- 
plays the strangeness content in the core and the crust 
of the neutron star as obtained in our calculations. The 
strangeness content is given by 

r s = xr? p + (l-x)rf P . (37) 

with 

r QP = )rf P = E^£ (38 ) 

where qf is the strange charge of baryon B. In both the 
cases, the strangeness fraction rises steadily. However in 
case of UQM the onset of strangeness starts at around 1.3 
Po and gets saturated at around 6.0 po where the matter 
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content in CFL starts at around 2 times normal nuclear 
matter density, then rises steadily in the mixed phase 
and ultimately gets saturated at ^6.6po- Using the EOS 
corresponding to the UQM and CFL we now compute 
the properties of hybrid neutron star similar to the case 
of pure neutron star mass, M Q and radius, gravitational 
red shift (Z), the baryonic mass of the star, Mb. 

In Fig. 5 the maximum mass of the hybrid star is 
plotted as a function of the radius of the star as obtained 
by the UQM and CFL model in our calculations. The 
radius of the maximum mass of the star is sensitive to 
the low density EOS. In order to calculate the radius 
and to plot it versus the star mass, we have used the re- 
sults of Baym, Pethick and Sutherland 45| for low bary- 
onic densities. The maximum mass obtained are 1.44 
M Q and 1.35 M Q for UQM and CFL respectively. The 
corresponding radius obtained are 9.35 km and 9.33 km 
respectively. These results are found to be in good agree- 
ment with different field theoretical models and also from 
observational point of view(4(| [4?J . The success of the 
model is evident from figure 7, where we plot M/R ratios 
as reported by [46ll47|. 

Constraints on the mass-to-radius ratio can be ob- 
tained from accurate measurements of the gravitational 
red shift of spectral lines produced in neutron star pho- 
tospheres. As already mentioned that a red shift of 0.35 
from three different transitions of the spectra of the X- 
ray binary EXO0748-676 was obtained [4y|, which corre- 
sponds to M/R — 0.15MQ/km. Another constraint to 
the mass-to-radius ratio given by M/R = 0.069M q //cto 
to M/R = O.H5MQ/km was determined from the obser- 
vation of two absorption features in the source spectrum 
of the IE 1207.4-5209 neutron star [47[. However, in the 
second case, the interpretation of the absorption features 
as atomic transition lines is controversial. The absorp- 
tion features are of cyclotron nature which make the re- 
lated constraints unrealistic [H, |49| . In our calculations 
the M/R ratio comes out to be 0.15 and 0.14 for UQM 
and CFL respectively. It is evident that our results are 
in very good agreement with the observed analysis. We 
have added the lines corresponding to those constraints 
shown by straight lines in Fig. 5. 

Baryonic mass as a function of maximum mass of the 
star is plotted in Fig. 6. It is to be noted that the baryon 
mass always exceeds the maximum mass, which is typical 
of the compact objects. The difference between the two 
is defined as the gravitational binding of the star. The 
baryonic mass, M},(Mq), obtained by UQM and CFL are 
1.63 Mq and 1.51 M Q respectively. 

Gravitational red shift of a star is an important prop- 
erty that is defined as 



TABLE II: Calculated properties of Neutron star 



Z 



1 



yjl - 2GM/Rc 2 



(39) 



EOS 


M(M ) 


E c {lO^gcm- 5 ) 


R(km) 


M b (M Q ) 


Z 


UQM 


1.44 


31.0 


9.35 


1.62 


0.35 


CFL 


1.35 


32.0 


9.33 


1.51 


0.32 



Red shift primarily depends on the M/R ratio of the star 
and observationally in absence of precise and direct mea- 
surements of mass and radius of the star it serves the 



7, we plot the gravitational red shift as a function of mass 
of the star with the straight lines as observational con- 
straints. The two solid straight lines are corresponds to 
red shift Z = 0.12 — 0.23 was determined from the ob- 
servation of two absorption features in the source spec- 
trum of the IE 1207.4-5209 neutron star [13] and the 
dotted dash line is correspond to Z = 0.35 which is de- 
termined from three different transitions of the spectra of 
the X-ray binary EXO0748-676 In our calculation, 
we obtained Z = 0.35 and Z = 0.32 for UQM and CFL 
respectively. Our results are in good agreement with the 
data. For the sake of completeness, we present the overall 
results in Table II. 



IV. SUMMARY AND CONCLUSIONS 

In the summary, we obtained the EOS of state for neu- 
tron star and studied the various properties of the star 
i.e., like mass, radius and redshift. The phase transitions 
from hadron to quark matter and the existence of the 
mixed phase in the inner core of the star are analyzed. 
The considered model based on the assumption that the 
neutron star consists of the inner core and crusts as two 
main parts of it. The theoretical model [13] for the outer 
region is described by the Lagrangian which includes the 
full octet of baryons interacting through the exchange 
of meson fields. The chosen parameter set based on ef- 
fective field theory motivated relativistic mean field in- 
cludes nonlinear scalar-vector and vector-vector interac- 
tion terms which leads to the soft EOS. This is in agree- 
ment of DBHF results. But the mixed phase EOS are not 
feasible, within the original G2 parameter set of the E- 
RMF formalism. To get a reasonable EOS for the mixed 
phase, for the hadronic matter, we changed the incom- 
pressibility marginally. 

The inner core of the star is described by UQM and 
CFL models. In UQM the quarks are treated as massless 
particles inside a bag of finite dimension and confine- 
ments results from the balance of the pressure on the 
bag walls from the outside and the bag pressure result- 
ing from the kinetic energy of the quarks inside the bag. 
The CFL phase consists of equal numbers of u,d and s 
quark and so requires no electrons to make it electrically 
neutral. In this paper we used different bag pressure 
B\= 170 MeV and 188 MeV for UQM and CFL model, 
respectively. 

The EOS are then employed to study and evaluate the 
global properties of the hybrid neutron star like mass, 
radius and gravitational redshift Z. The masses are pre- 



7 



9.33 km in E-RMF+UQM and E-RMF+CFL models, re- 
spectively. We have also shown that the EOS with UQM 
and CFL quark phase satisfies the constraint imposed 
by the recently measured redshift of 0.35 from three 
different transitions of the spectra of the X-ray binary 
EXO0748-676 0. The redshifts are 0.35 and 0.32 in 
the E-RMF+UQM and E-RMF+CFL formalism, respec- 
tively, which coincides very well with the measurement. 
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FIG. 1: Equation of state obtained with E-RMF model plus UQM (solid line) and plus CFL (dash-dotted line). 
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FIG. 6: Baryonic mass as a function of maximum mass of the neutron star with E-RMF plus UQM (solid line), E-RMF plus 
CFL (dashed line). 
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FIG. 7: Gravitational Redshift (Z) as a function of maximum mass of the neutron star with E-RMF plus UQM (dashed line), 
E-RMF plus CFL (dotted line) with experimental observations (solid and dash-dotted lines) . 



